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As part of a larger project to estimate the fs decay constant, we are recalculating f^ attc using a variational 
smearing method in an effort to improve accuracy. Preliminary results for the static Bb parameter and HQET 
two point functions are also presented. 



1. INTRODUCTION 

The extraction of CKM matrix elements from 
experimental data requires the calculation of 
QCD matrix elements of operators involving the 
bottom quark jlj]. For a number of years the 
MILC collaboration has been doing a systematic 
study of the Jb decay constant Q . The difficulty 
of simulating the bottom quark, with its mass 
larger than current feasible inverse lattice spac- 
ings, was overcome by interpolating between the 
results of simulations using Wilson quarks with 
masses around the charm mass and those from 
static simulations (infinite mass). 

Unfortunately, for some of the simulations the 
jstatic resu n s were no t useable. These results 
came as a by-product of the hopping parameter 
expansion of the heavy quark, and for technical 
reasons had significant contamination by higher 
momentum intermediate states when the physical 
volume was large J2| . In addition, it is well known 
that the poor signal to noise ratio of the static 



simulations makes it important to use an efficient 
smearing method. To provide f^ atlc results on 
all lattices, and to reduce the errors even in the 
cases where the results were previously available, 
we have started a set of static-light simulations 
using a variational smearing method. 

2. STATIC f B 

The computational cost of our current imple- 
mentation of the FFT on parallel machines, pro- 
hibited the use of a more sophisticated smearing 
technique such as MOST so we use the fol- 
lowing basis of smearing functions: 



s(r)i = e- A *\^\ 
P ~M r I 



(l-B- 



♦Presented by C. McNeile. 



s(r) 2 = ( 
s(r) 3 = e~ A *l - U(l-C*|r|--D*|r| 2 ) 

The parameters A, B, C and D were obtained 
from uncorrelated fits to the Kentucky group's 
measured wave functions of static-light mesons 
(at = 6.0) f|. The A, B, C and D param- 
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Table 1 



Simulation parameters 





volume 


# configs 




#« 


5.6 


16 3 x 32 


100 


0.01 


3 


5.5 


24 3 x 64 


100 


0.1 


3 


5.445 


16 3 x 48 


100 


0.025 


3 



eters were scaled to the appropriate value for a 
given j3, using the estimates of the lattice spac- 
ing. Because of uncertainties in the lattice spac- 
ing, and to have the flexibility to choose different 
smearing functions for each kappa value, two ad- 
ditional sets of the parameters were chosen. Thus 
a variational smearing matrix of order ten (in- 
cluding the local operator) is used in the simula- 
tions. The static quark is smeared relative to the 
light quark using standard FFT methods. The 
completed static-light production runs are shown 
in Table [l]. The (3 = 5.6 configurations were gen- 
erated by the HEMGCC collaboration f|. Our 
analysis of the data is very preliminary-in partic- 
ular we have not yet fully optimized our smear- 
ing functions. All the results presented here use a 
single exponential source. We will focus on com- 
paring the new static results with the numbers 
from other simulations. Ultimately, all the static 
data will be combined with that from propagating 
quark simulations [Q. 

We do a simultaneous fit to the smeared-local 
and smeared-smeared correlators. To compare 
our raw lattice numbers with other simulations, 
we quote our numbers in terms of Zl, defined by 

C LS (t) = Z L Z s e- E "^ (1) 
C SS (t) = ZsZse-^ (2) 

Zl is related to the decay constant, no pertur- 
bative factors are included in its definition and 
we assume the light quark propagator has been 
multiplied by 2k. E S i m is the energy of the static- 
light meson, it is equal to the sum of the difference 
between the mass of the B meson and the bottom 
quark mass, and an unphysical - renormalization 
factor. Correlations were included for the fits in 
time, but no kappa correlations were included in 
the chiral extrapolations. 



Table 2 

Zl fit results for (3 — 5.6, am sea = 0.01, with a 
source exp (—0.4 * | r |), and fit region 7 to 11 



K 




a 3 ' 2 Z L 


X 2 /dof 


0.156 


0.708^ 


0.312±4 


4.7/7 


0.158 


0.684^ 


0.287±g 


4.9/7 


0.159 


0.637^4 


0.274^ 


5.6/7 


0.16103 


0.649±| 


0.249i£ 


0.01/1 



Some preliminary, static fit results for the = 
5.6 simulation are contained in Table ^. Ali Khan 
et al. || have also calculated Zl on gauge con- 
figurations from the same (3 — 5.6 simulation, 
as those used here, but not on exactly the same 
sample of configurations. At K = 0.1585, they 
quote aE slm = 0.528(5) and a 3 / 2 Z L = 0.24(3). 
Although we do not have this kappa value in our 
simulation we can use the information obtained in 
the chiral fit model to estimate a 3 / 2 Z L = 0.280tg 
and aE S i m = O.536I4, (where we have added 
loguo, with uq = 0.867, to our value of aE S i m 
because Ali Khan et al. rescaled all their gauge 
fields by the tadpole improvement factor of uo). 
The results from the older MILC static calcula- 
tion also agree with those from the new simula- 
tions. 

We have done some fits to the (3 = 5.445 static 
data. The masses obtained were consistent with 
the older MILC calculation. This is a large- 
volume case where the older method does not 
produce useable static Zl factors, so no check 
is available there. However, the new static corre- 
lators seem reasonably consistent with the prop- 
agating quark results. 

3. STATIC B B PARAMETER 

The Bb parameter is required in the extraction 
of the Vtd CKM matrix element from the experi- 
mental data on B—B mixing. For the (3 — 5.5 and 
5.6 simulations we calculated the static Bb pa- 
rameter (so far we have only analyzed the (3 = 5.6 
data). Ours is the first calculation of the static 
Bb parameter that includes dynamical fermions. 
The method used is described in references B 
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Figure 1. Bl parameter as a function of time 
slice, for k — 0.156, and source exp (—0.4 * | r |) 

and . The same set of smearing functions used 
in the fp atic calculations was used to smear the 
quarks in the external mesons. Fig. [l] shows the 
static Bl operator as a function of time. In Ta- 
bic H we show some preliminary results for static 
Bb parameter. The x 2 /dof for all the fits in Ta- 
bic H were all close to 0.5. The errors are statisti- 
cal only; the systematic errors due to the choice 
of fit range are larger than the statistical errors. 
The errors should be reduced when we include the 
additional smearing functions in our analysis. 

We used the Kentucky group's Q organization 
of perturbation theory to find the required linear 
combination of operators to calculate 
However we omitted some next to leading or- 
der log/x/m terms. The missing terms have only 
recently been calculated B and are a small ef- 
fect. Chirally extrapolating the Bsims) results 
to k c — 0.16103 and converting the results to the 
one loop RG invariant Bb parameter, we obtain 
the preliminary value of Bb = 1-31+| (statisti- 
cal error only). This result is consistent with the 
values obtained from quenched simulations using 
Wilson fermions ||. 

4. LATTICE HQET 

To calculate the form factors of the semi- 
leptonic decays of the B meson, we want to follow 
a strategy similar to the one used in the fs simu- 



Table 3 

Static Bb parameter results for (3 = 5.6, am sea = 
0.01, using the source exp (— 0.4 * | r |), fixing 
timeslice 29, and fitting times 2 to 6. 



0.156 



0.158 



0.159 



B L 


0.98(1) 


0.99(2) 


1.00(2) 


Br 


0.96(1) 


0.94(2) 


0.94(2) 


Bn 


1.01(2) 


1.01(2) 


1.01(2) 


-IBs 


1.00(1) 


1.00(1) 


1.01(2) 


B B (m B ) 


0.91(2) 


0.92(2) 


0.93(2) 



lations, except that we will combine the results of 
heavy quark effective field theory (HQET) with 
the analogous propagating quark calculations, to 
allow the final results to be interpolated to the B 
meson mass. 

As a "warm up exercise" we studied the two 
point function of a HQET-light meson. This ex- 
ercise allows us to investigate the smearing of the 
quarks in the B meson and to study the nonper- 
turbative renormalization of the velocity (a pecu- 
liarity of lattice HQET) - both of which are im- 
portant prerequisites to the calculation of form 
factors. 

We have implemented the HQET propagator 
equation introduced by Mandula and Ogilvie || . 
The two-point function for an HQET-light meson 
at finite residual momentum p and bare velocity 
v is 

C(p,t;v) = J2J2f^y PX 

X r 

(K(r, 0)759(Q, 0)q{x, t) l5 b v (x, i)> 

When there is no excited state contamination, the 
correlator has the form 



C(p,t;v) = Z v f Z v L e- E{p ' vR)t 



(3) 



where ZJ and Z v h are the smeared and local ma- 
trix elements. The dispersion relation for an 
HQET-light meson is M] 



E(p,v R ) = 



v R .r> 



(4) 



where E sim is equal to the energy obtained in the 
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Figure 2. Effective mass piot for the HQET- light 
two point function, with source exp (—0.67 * | r |) 

static simulation. Here v R is the renormalized 
velocity. 

As a pilot study we generated 20 correlators 
at P = 5.445 with a light quark kappa value of 
0.160. So far, a single exponential source, as in 
the static simulations, was used. Fig. || shows 
the effective mass plots for a HQET-light meson 
at zero residual momentum with bare velocities 
only in the x direction of 0.1, 0.5, and 0.8; the 
static effective mass plot is also included (equiv- 
alent to zero velocity). The single exponential 
produces almost usable plateaus for all three ve- 
locities - the plateaus should improve when we 
use variational smearing. The signal to noise ratio 
does not decrease rapidly with increasing veloc- 
ity. The renormalization of the velocity is caused 
by the breaking of the Lorentz symmetry by the 
lattice. This renormalization is required in the 
calculation of form factors such as the Isgur-Wise 
function |0. Our preferred way to extract the 
renormalization is to use the dispersion relation 
in Eq (Q) with the energy taken from simulations. 
This approach seems to us |l^] Q to be simpler 
than the one used by Mandula and Ogilvie [|l3| . 

An estimate of the velocity renormalization can 
be obtained from the data in Fig. |^, using naive 
fits to the effective masses and using Eq Q with 



zero residual momentum. For the bare velocity 
of 0.5 (0.8) in the x direction, the renormalized 
velocity was approximately eighty (sixty five) per- 
cent of the bare velocity. With higher statistics 
and better smearing a more sophisticated analy- 
sis will be done. The perturbative calculation of 
Mandula and Ogilvie jl3| also gave a renormal- 
ized velocity smaller than the bare velocity. 

This work is supported in part by the U.S. De- 
partment of Energy and the NSF We would like to 
thank Terry Draper for discussions on HQET and 
Bb, and Joe Christensen for providing us with 
the coefficients for the Bb parameter calculation. 
The runs are being done on the Cornell Theory 
Center's SP2, and on the 512 node Paragon at the 
Oak Ridge Center for Computational Science. 
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